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We make a detailed eikonal fit to current data on the total and elastic scattering cross sections, 
the ratios p of the real to the imaginary parts of the forward elastic scattering amplitudes, and 
the logarithmic slopes B of the differential cross sections da/dt at t = 0, for proton-proton and 
antiproton-proton scattering at center-of-mass energies W from 5 GeV to 57 TeV. The fit allows 
us to investigate the structure of the eikonal amplitudes in detail, including the impact-parameter 
structure of the energy-independent edge in the scattering amplitude shown to exist by Block et al. 

[I[. We show that the edge region has an essentially fixed shape with a peak at approximately the 
“black disk” radius Rtot = \Jatot iT-it of the scattering amplitude, a constant width tedge ~ 1 fm, and 
migrates to larger impact parameters with increasing energy proportionally to -Rtot- We comment 
on possible physical mechanisms which could lead to the edge. We show that the eikonal results for 
the cross sections and p values are described to high accuracy by analytic expressions of the forms 
used in earlier analyses by Block and Halzen, and extend the result to the elastic-scattering slope 
parameter B. These expressions provide simple extrapolations of the results to much higher energies 
where the cross sections approach the black disk limit with Ueias, Uinei —>■ crtot/2 and B —>■ (Ttot/87r. 
Finally, we calculate the survival probabilities for large rapidity gaps in the scattering. 

PACS numbers: 13.85.Dz,13.85.Lg, 13.85Tp 


I. INTRODUCTION 

In a recent paper [ij Block et al. established that the proton-proton scattering amplitude has an edge in impact 
parameter space with a width that remains essentially constant over many orders of magnitude in the center-of-mass 
energy W. This result was derived using the general forms of the scattering amplitudes in impact parameter space for 
strongly absorptive scattering, and the very accurate Block-Halzen ht to the proton-proton (pp) and antiproton (pp) 
total and inelastic cross sections and ratios p of the real to imaginary parts of the forward elastic scattering amplitudes 
for 6 < IF < 1800 GeV Q. That ht incorporated the asymptotic In^ s limit of the growth of the cross sections at 
large s = W^, the constraints on the phase of the scattering amplitude imposed by analyticity and crossing symmetry, 
and constraints on its magnitude and slope at 4 GeV implied by consistency with low-energy data. It successfully 
predicted the cross sections subsequently measured in the multi-TeV range at the Large Hadron Collider (LHC) and 
in cosmic ray experiments ii- 

The results in [I| did not depend on a particular ht to the scattering amplitudes, but only on the general forms of 
the amplitudes in impact parameter space. In the present paper, we use a detailed eikonal description of the scattering 
to ht the pp and pp data on dtot, CTeias, P, and on B, the logarithmic derivative of da/dt at t = 0. The ht includes 
a number of new measurements at higher energies. While the results we obtain for these quantities are essentially 
equivalent to those obtained earlier, the detailed ht allows us the investigate the structure of the eikonal amplitudes 
including, in particular, the structure of the edge in impact parameter space. 
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In the following sections, we first establish our conventions and give expressions for the cross sections, p, and B in 
terms of the eikonal function (Sec.|ITl). We use a general parametrization of the eikonal function which, importantly, 
incorporates the power-law growth oc s'^ found in QCD-based minijet Q and Reggeon i0 models, and the exponential 
cutoff in impact parameter suggested by the proton form factor. This results in the asymptotic approach of the 
scattering amplitudes to the black-disk limit in which utot c)c In^ s, Ueias, crinei —t fftot/S, B —>■ utot/STr, and p —?► 0. 

We present the results of our fit to the pp and pp data in Sec.lIIIl We then use the results to investigate the structure of 
the eikonal amplitudes ISec. lIV XI) . including the relative importance of different contributions to the eikonal function 
and the (slow) approach of the scattering to asymptotic behavior dominated by gluon-related processes. 

We show in Sec. lIVBi that the description of the total and elastic cross sections and p values obtained in the eikonal 
model can be fitted to high accuracy by analytic expressions of the form used in the Block-Halzen fits to the data 
justifying their assumptions. The resulting expressions for these quantities can be extrapolated reliably to 
higher energies. We also extend this analysis to B, and give results that may be useful in other contexts such as the 
analysis of cosmic ray cross sections. We prefer these extrapolations to those using the eikonal model for reasons we 
discuss. 

We then investigate the structure of the edge region in the scattering amplitude ISec. lIV Cll . We find that the edge 
maintains a nearly constant shape in impact parameter space, with a width which remains essentially constant at ^ 1 
fm up to the highest energies studied to date, and presumably to much higher asymptotic energies. We comment on 
some possible explanations of the edge and its form at large impact parameters in Sec. lIVE)! this is a problem that 
needs further study. 

We note in Sec. lIVE)] that the component of the cross section associated with the edge gives the Pumplin bound 
Q on single-particle diffraction dissociation in pp or pp collisions, and suggest that experiments to test the bound 
would be useful. Finally, in Sec. lIVFl we discuss and calculate the survival probabilities for large rapidity gaps in the 
scattering. The details of our eikonal model are discussed in the Appendix. 


II. CONVENTIONS 

In the following, we will be concerned with proton-proton (pp) and proton-antiproton (pp) scattering at high 
energies. We will neglect the (presumably small) effects of the nucleon spins, and describe the scattering amplitude 
and cross sections in an impact parameter or eikonal representation; this is valid at small angles when many partial 
waves contribute to the scattering and the (unitary) partial wave series can be converted to an integral over the 
impact parameter. We will write the resulting spin-independent eikonal scattering amplitude and differential elastic 
scattering amplitude as 


f{s,t) = Jo(6\/=t), (1) 

^(s,t) = 7r|/(s,t)|^ (2) 

Here s = = 4(p^+m^) is the square of the total energy in the center of mass (c.m.) system, p is the c.m. momentum 

of either incident particle, b = j/p where j is the partial-wave angular momentum, and t = —2p^(I — cosd) is the 
invariant 4-momentum transfer for elastic scattering at the angle 9. We will define the eikonal function x(6, s) as 
X = XR + note that some other papers use different conventions, e.g., nil- 
With these conventions, the elastic, total, and inelastic cross sections are 
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The ratio p of the real to the imaginary part of the forward scattering amplitude and the logarithmic derivative B 
of the differential elastic scattering cross section at t = 0 are also frequently measured and will be used in our analysis. 
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An accurate approximation for B when the real part of the scattering amplitude is small is to set XR = 0- Then 
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We have used the exact expression in Eq. ([SI) in fitting the experimental data, but note that the approximate expression 
would have been adequate. 


III. FIT TO HIGH ENERGY PROTON-PROTON AND ANTIPROTON-PROTON SCATTERING DATA 


The model we have used in fitting the high evergy pp and pp cross sections is a modification of the “Aspen model” 
of Block et al. HP which was motivated by the structure of the eikonal function found in QCD minijet models 
for the scattering. We will follow the notation used in those references even though the precise identification of the 
terms made there cannot really be maintained in a more general setting. We write the eikonal functions in terms of 
crossing-even and crossing-odd components, with 


Xpp{b,W) = [xE{b,W) + xo{b,W)]/2, 
Xpp{b,W) = [xE{b,W) - xo{b,W)]/2. 


The even and odd functions are defined as 


XEib,W) = i 
Xo{b,W) = -CsEgg 


2—20:1 
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where the phases of the functions in Eqs. m and m are determined by the constraints imposed by analyticity and 
crossing symmetry [§, [ll| . 

In these expressions, the factors A{b, p) are overlap functions for the colliding hadrons and the “cross sections” 
CTy are intended to describe the interactions between the the corresponding components i and j of the two particles 
chosen from the matter (g) or gluon (g) fields. The details of the model are given in Appendix lAl 

Our parametrization of x is general and very flexible, including a leading power-law dependence additional 
logarithmic and constant terms, and falling Regge-like terms in s. Our objective is to get a good fit to all the pp 
and pp data up to the highest energies where measurement exist, and to then use the results to study the eikonal 
structure of the scattering amplitudes with immediate emphasis on the edge region [l|. In contrast to our relatively 
free parametrization of y, other recent parametrizations such as those in [5l-l7| are based on specific dynamical models, 
and those papers emphasize the testing of those models through fits to the data. 

We have used our parametrization and the expressions above to fit the combined data on pp and pp total cross 
sections for W > 5.3 GeV and the elastic scattering cross sections, p, and B for energies W > 10 GeV. The fit was 
further constrained as described in Q by fixing the values of (itot.pp and (Ttot,pp at W = 4 GeV to match the results 
obtained from the extensive low-energy data. This is the same energy range with the same constraints as used in 
the Block-Halzen fits based on analytic amplitudes with a In^ s high energy behavior ii- However, we include the 
newer data at very high energies from the Large Hadron Gollider (LHG) |12l - [l^ and the Auger [13 and HiRes (l^ 
collaborations. 

The fits were performed using the sieve algorithm [l^ to eliminate 13 outlying points among 179 total datum points. 
Nine parameters were used in the fit leaving 157 degees of freedom, a total x^ of 173.0, and a raw y^/d.o.f. = 1.10. 
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FIG. 1: Top panel: fits to ntot.pp (blue dots and solid line) and ntot.pp (red squares and dashed line). Only data above 5 GeV 
were used in the final fit, with the cross sections constrained to fit compilations of low-energy data at 4 GeV [^. . Bottom 
panel: fits to (Teias.pp (blue dots and solid line) and (Teias,pp (red squares and dashed line). The fit used only data above 10 GeV. 


This must be renormalized by the sieve factor TZ k, 1.1 to 7?.x^/d.o.f. = 1.21 to account for the elimination of the 
outliers M- The total would increase by 113.6 if we included the outliers, so the change would be substantial. 
For comparison, the x^/d.o.f. given by the fit is just 1.15 for the pp and pp total cross sections and p values alone; 
much of the increase in the final result comes from the fit to the rather scattered values of B. We note that all datum 
points including the outliers omitted in the final fit are shown in the figures comparing the fits with data. 

The results for the fits to the total and elastic scattering cross sections are shown in Fig.[T] The fits to the p values 
and the logarithmic slopes B of the forward differential elastic scattering cross sections da/dt, Eq. ([7]), are shown in 
Fig. [2] The highest energy data for p are from the LHC at 1,800 GeV. The value predicted for the LHC at IF = 7 
TeV is p = 0.133. The data for B include the TOTEM results [H, [3 from the LHC at VF = 7 TeV. 

The measured and predicted differential cross sections da/dt are shown in Fig. [3] at IF = 1800 and 7000 GeV. 
Our descriptions of the cross sections at small |t| are good, corresponding to our fits to the B parameters, and the 
locations of the diffraction minima are reproduced properly. We are not concerned about the failure of our simple 
eikonal model to reproduce the differential cross sections in detail at large values of |t| since the scattering amplitudes 
in this region are very sensitive to the cancellations which result from the oscillations of the Bessel function in Eq. o, 
with the resulting scattering amplitudes of order ^ 10 ^ x /(s,0). As an illustration, we show the integrands for 
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FIG. 2: Top panel: fits to the ratios p of the real to the imaginary parts of the forward scattering amplitndes for pp (blue dots 
and solid line) and pp (red squares and dashed line) scattering. The horizontal dashed line is at p = 0. Bottom panel: fits 
to the logarithmic slopes of the elastic differential scattering cross sections da/dt for pp (blue dots and solid line) and pp (red 
squares and dashed line) scattering. 

Im/(s,i) for ^/s = W = 1 TeV and |t| = 0.5 GeV^ and 1 GeV^ in Fig. SI The existence of large cancellations and 
the resulting sensitivity of the integrals to small details of the eikonal function not modeled here are evident. We 
emphasize, however, that the cross sections, p, and B are much less sensitive to such details. 


IV. STRUCTURE OF THE EIKONAL AMPLITUDES AND THE EDGE IN THE pp AND pp 

SCATTERING AMPLITUDES 

A. Eikonal structure 

It will be important for later interpretation to understand the relative importance of the various contributions to 
the eikonal function and cross sections. Since XR is small, xi determines the cross sections to good approximation. 
In the top panel of Fig. [5] we therefore compare the imaginary parts of the energy-dependent factors in xi- The most 
important contribution at high energies (e.p., W > 1 TeV) is clearly that from gluon-gluon scattering, parametrized 
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FIG. 3: Top: the differential cross section da/dt from the E710 experiment [^. at VF = 1800 GeV. Bottom da/dt from the 
TOTEM experiment at IF = 7000 GeV. 


in our model in a form suggested by the rapid growth of gg scattering in perturbative QCD and the corresponding 
minijet models for the rise in dtot with energy. 

The other important contribution at high energies is that labelled Uqq. This term cannot be separated in the fit 
from Gqg, and is not to be interpreted strictly in terms of quark-quark scattering in the sense of the parton model. 
The qq scattering in that model in fact becomes small at lower energies, while the combined contributions of Uqq and 
aqg increase as parametrized. The low-energy behavior presumably arises from “soft” processes such as the scattering 
of valence quarks evident, for example, in the approximate 2/3 ratio of the irp and pp cross sections, and to Regge 
exchange terms. 

The crossing-odd contribution xo vanishes rapidly with increasing energy. Finally, the mixed quark-gluon term, 
parametrized in a form suggested by minijet models [lOl| , is strongly mixed and correlated with the other terms in the 
fitting, and should not be interpreted directly in terms of qg scattering. 

The eikonal factors and the complete integrands in the expressions for atot, Cinei, and Ceias in Eqs. GD-® are shown 
for pp scattering at a progression of energies in Fig. [SI The quantities rj and cr used in the labels in this figure are 
p = and cr = cosxr] similarly, sr = siny/i. 

We note several important features of the curves shown. First, the eikonal factors {1—CRp) for ctot and {l—2cRp+p^) 
for cTeias he well below the asymptotic “black disk” limit 1 at small values of the impact parameter b for energies 





7 



FIG. 4: Plot of the integrand 6(1 — cosxne ^^)Jo(fe%/^) for the imaginary part of f{s,t) versus the impact parameter b for 
W = 1 TeV and \t\ = 0.5 GeV^ (solid blue curve) and 1 GeV^ (dashed purple curve). 


W < 5 TeV as shown in the left-hand column in Fig. [51 The scattering is far from asymptotic, and the approach to an 
asymptotic distribution flat at the value 1 out to a sharp cutoff radius ~ i? is extremely slow. The inelastic integrand 
approaches asymptotic-like behavior more rapidly, with the eikonal factor 1 — tv 1 becoming flat at 1 at small 6, at 
a sharper cutoff, at lower energies. 

The actual integrands including the geometric factor b are shown in the right-hand column in Fig.|6l This factor 
pushes the relevant impact parameters toward larger b and introduces the peaked behavior shown. The main integrals 
involved in the calculation of the logarithmic slope parameter B = (i[ln {da/dt)]/dt\t=o, Eq. ([8]) or ([9|), involve an extra 
factor b^ in the numerator, with the result that the main contributions to the integral are pushed to larger values of 
b and become increasingly sensitive to the tail of the eikonal distribution for deias- 

The integrand for p, the ratio of the real to the imaginary parts of the forward elastic scattering amplitude, Eq. ([6|), 
involves a factor rj = siny/{ in the numerator. Since the transparency factor rj vanishes strongly at small b at 
high energies, the result of strong inelastic absorption, the main contributions to p are pushed toward higher values 
of b, beyond the peak in the integrand for (Xeias, and are again more sensitive to the tail of the distribution than the 
elastic scattering cross section itself. In addition, the phase of the scattering amplitude, hence the ratio of XR to %/, 
is determined at high enercies by the form assumed for x coupled with the constraints imposed by the analyticity of 
the scattering amplitude iini , xr is not freely variable. As a result, there is a tension between B and p when fitting 
data: both are sensitive to the tail of the distribution and the two parameters are therefore coupled with respect to 
changes in that distribution. The fit to the current data discussed above incorporates the constraints imposed by 
both B and p. 


B. Connection to the description in terms of real analytic ampitudes 

Despite the seeming complication of our fit to the data in terms of the eikonal function, the results we obtain for 
CTtot and p for pp and pp scattering can be described very well with an expression of the simple form used by Block and 
Halzen [2, Q in their earlier ht to the corresponding data up to FF = 1800 GeV. That fit gave successful predictions 
of the more recent, higher energy data ii- 

The Block-Halzen analysis assumed a In^ s bound on the growth of the cross sections at high energy, imposed the 
constraints implied by the analyticity of the scattering amplitudes, and was constrained to connect smoothly to the 
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W, GeV 

FIG. 5: Top: Comparison of the imaginary parts of the different energy-dependent factors in the eikonal function. Solid (red) 
curve: o-,,. Dot-dashed (blue) curve barely visible near zero amplitude: cr^g. Long-dashed (purple) curve: Ogg. Dotted (black) 
curve: the odd term. Bottom: comparison of the cross sections calculated with (solid blue curve) and without ( dashed red 
curve) the inclusion of the gluon-gluon {gg) term in the eikonal function. The cross section for pure gluon scattering is shown 
as the long dashed purple curve. 


low-energy data. It was based on the use of analytic amplitudes of the form 


--(y) = co+c,\n(-)+C2ln^(-)+P'(-y ' ± S (-) 
\mJ \m/ \m/ \m/ 


(7=*= 


{|ci+7rc2ln(^)-/?'cot(^)( 


47r /TTQ: 

+ -/+(0)±<5tan(— 




a-1 


(14) 

(15) 


where the upper and lower signs are for pp and pp scattering respectively. Here v is the laboratory energy of the 
incident particle, with 2mv = s — 2m? = — 2m? where m is the proton mass. 

Their fit used the then-extant data on ctot and p in the range 6 < VF < 1800 GeV plus analyticity constraints on 
the values and slopes of the cross sections at TV = 4 GeV which followed from finite-energy sum rules applied to the 
data at lower energies. The fit was excellent and led to successful predictions of the results for the cross sections later 
measured at the LHC and in cosmic ray experiments ii- A later analysis of the inelastic scattering data using an 
expression of the same form as Eq. (fl3 was also successful and gave evidence of an approach to the black disk limit 
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FIG. 6: Plots of the eikonal factors (left-hand column) and those factors multiplied by the geometric factor b in the integrands 
for the pp total cross section (top row), inelastic cross section (middle row), and elastic cross section (bottom row). In the 
labels for the ordinate, rj = and ca = cosxfl- The curves in each panel correspond, bottom to top, to energies IT = 50 
GeV (red curve), 500 GeV (brown curve), 5 TeV (blue curve), 50 TeV (purple curve) and 1000 TeV (black curve). 


at ultrahigh energies i , with (Tinel -t Crtot/2. 

We have checked that the use of the Block-Halzen expressions to fit “data” derived from our results gives curves 
for tTtot and Ueias that are almost indistinguishable from the curves in Fig.[TJ Both fits describe the data quite well, 
and we conclude that they are consistent. 

The expressions in Eq. (1141) and Eq. (TTSl) simplify to a more familiar form for m?, with vjm ^ j2m? = 

sl2rr?. The corrections to the logarithmic terms are negligibly small for W in the region of the fit. The corrections 
to the power-law terms are a fraction of a millibarn for FT = 4 GeV, negligible for FT > 6 GeV, and can be absorbed 
overall in slight adjustments of the powers and coefficients of those terms where they are relevant. As a result, the 
argument v in the formulas for (7^(i^) and p^{v) can be converted directly to sl2rr?^ or with some rearrangement of 
terms and coefficients, to Wjm^ without loss of accuracy in the region used in our fit and that of Block and Halzen. 
A similar expansion quadratic in ln(s/m^) for CTeias follows from the results in Q. The coefficient of In^ (s/m^) in CTinei 
was found in Q to be one-half that found for utot as required for an asymptotic black-disk limit for the scattering 
amplitude. 

A quadratic in Ins was also used by Schegelsky and Ryskin (2^ to fit the data on B. Their result for the coefficient 
of In^ s was consistent with the Block-Halzen fit to the cross sections and p alone and the expectation for black disk 
scattering that B Rtot/4 = CTtot/87r for s —>■ oo. 

Similarly, our results for Bpp can be written to an accuracy of a few parts in 1000 in the same form as 

Bpp{W) = 7.229-b 1.0862 ln(FT/mo)-b 0.02209ln^(FT/mo)+3.719(mg/FT2)i/2 GeV"^ (16) 

for 6 GeV < FT < 10^° GeV, where mg = 0.6 GeV is our scale factor and FT = -^s. This result is potentially useful 
in the conversion of cosmic ray cross sections for proton-air scattering to pp cross sections. 
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While our eikonal fit to the complete data set gives results for the cross sections and p that are essentially equivalent 
numerically to those of Block and Halzen over the energy range currently accessible, it is not immediately clear 
analytically from the rather complicated eikonal expressions why the simple expressions in Eqs. 0 and m — or 
their reduced high-energy forms — should work so well. However, the asymptotic In^ s growth of the cross sections 
in Eq. (03 and the approach to the black-disk limit of the scattering follow directly from the expected power-law 
growth of the eikonal function with s coupled with the exponential cutoff in the overlap functions A{b, A), Eq. (jASp . 
for Xb ^ 1, an argument familiar in discussions of the Froissart bound. The first leads to strong growth of x/(6, W) 
with y/ ^ 1 and <C 1 at large s and small 6, and a corresponding saturation of the scattering amplitudes for 
CTtot, CTinei, and CTeias at the value 1 as seen in the left-hand column in Fig. [6] This saturation persists out to values 
of b such the exponential decrease in A{h, A) pushes xi to values below 1 beyond which the scattering amplitudes 
vanish exponentially. The rough condition X/(6, VE) < 1/2 determines the effective cutoff radius in b which, given the 
exponential behavior of A{b, A), can grow only as Ins. 

Simple arguments using Eqs. dS])-® and dH) or s then show that Ueias and tTinei tend to CTtot/2 for s sufficiently 
large, while B —^ ntot/STr, with all proportional to In s plus logarithmic and lower-order corrections associated with 
the edge region in the amplitudes, the region around the peaks in the integrands shown in the right-hand column of 
Fig.m These arguments provide a justification for the Block-Halzen form for the cross sections at sufficiently high 
energies; detailed checks using our eikonal fit show that the simple quadratic expressions are accurate at present-day 
energies. Finally, one can show from Eq. dH) and the constraint on the phase of the scattering amplitude imposed by 
analyticity and unitarity iEl that p —>■ 0 at very high energies. 

We emphasize that the eikonal fit allows us to calculate important quantities such as B and the differential scattering 
cross sections da/dt that are not accessible through a Block-Halzen type analysis without further input. Our analysis 
of the eikonal results on B shows that it, like the cross sections, can be described to high accuracy by a quadratic in 
Ins plus low-energy Regge-like terms, thus providing the necessary input. Given the uncertainty in the rate of the 
power-law growth of the eikonal function, and the uncertainty in its form at large b discussed in Sec. IIVDI we believe 
that fits to the data using the quasi-universal high-energy expressions in Eqs. da and da, the corresponding result 
for tTinei in i, and the expansion above for B, are likely to give a more reliable way at this point of extrapolating the 
cross sections to ultrahigh energies . 


C. The edge of the pp and pp scattering amplitudes 

Block et al. [l| recently established that the proton-proton scattering amplitude in impact parameter space has an 
edge region the width of which is essentially constant over many orders of magnitude in the center-of-mass energy 
W. This result followed from the usual form for the scattering amplitude combined with the assumption that the 
scattering is strongly absorptive. In particular, it was shown in [l| that this edge could be isolated using the properties 
of the transparency function r]{b, s) = exp [—%/(&, s)]. This function is very small at small impact parameters where 
the absorption is strong, then rises to unity — complete transparency and no scattering — at large impact parameters. 

This observation was exploited in [ij by noting that 


f^tot 2(7elas — 2(Tinel ^tot 

(17) 

= 47r / db br} {cos xr — v) 

Jo 

(18) 

poo 

dbbr]{l — if) 

Jo 

(19) 


for small real parts of the scattering amplitude, a condition satisfied in the present fit. The integrands in these 
expressions have the property that they are large only in the transition region between strong absorption and no 
scattering. 

Since Utot , CTeias, and CTinei are measured quantities, experiment gives a direct measurement of the edge integral in 
Ea J|17|) . To obtain its extrapolation to very high energies. Block et al. [I| used the very accurate Block-Halzen fit 
d, d to the pp and pp total cross sections and p values for 1800 > W > 6 GeV, and its extension to (Tinei This 
fit, which incorporated the asymptotic In^ s limit on the growth of the total cross sections for large s = and the 
constraints imposed by the analyticity of the scattering amplitudes and the lower energy data, successfully predicted 
the recent LHG and cosmic ray results Si- The constancy of the edge width followed directly from the use of those 
results in Eq. (II3, and did not depend on the detailed impact parameter distribution in Eq. (1191) . 

In the present eikonal fit to the pp and pp data, the edge integrand br]{cosxR — v) ^ br]{l — rj) is peaked at values 
of the impact parameter somewhat beyond the peak in the integrand for utot as shown in Fig. [3 and well into the tail 
region in the eikonal distribution for ctot as can be seen by a comparison to Fig.[6l top left. Not surprisingly, this is 
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just the region that determines the effective black disk radius Rtot = \/crtot/27r of the scattering amplitude. As seen 
in the comparison of the actual pp scattering amplitude with the black disk amplitude with the same value of utot in 
Fig.0 the “missing” contributions to the black disk amplitude for b < i?tot are supplied by the tail of the distribution 
with b > Rtot, with i?tot corresponding very closely to the peak in the edge integrand. 



FIG. 7: Comparisons of the integrands b(l — cosyrse for crtot (solid red curve), the “black disk” integrands for disk radius 
Rtot = -^(Ttot/STT (long-dashed black curve), and the edge integrands 6(cosyH — (short-dashed blue curves) at 

energies W — I, 5, and 50 TeV, left to right. The vertical and horizontal scales give the integrand and b in fm. 

The value of the edge integral should be approximately the height of the peaked integrand times its width tedge at 
half maximum. We define tedge as the edge width. Since 77(1 — 77) has a maximum value of 1/4 and b is approximately 
equal to Rtot at the peak, 


^tot 2fJelas ~ '^Rtottedge, (^tl) 

or 

tedge ~ (^tot 2(Telas)/'\/(^t) 

The edge width tedge was evaluated in [l| using the Block-Halzen fit to the pp and pp total cross sections and p 
values @ and its extension to the inelastic cross sections Q. It was found to be remarkably constant at tedge ~ 1 
fm above about 10 GeV. The result obtained here using our fit to the pp and pp data is essentially the same; this as 
shown in Fig.[8l 



FIG. 8: Plots of the pp edge width tedge calculated using the present eikonal fit to the pp and pp data (solid blue curve) and 
the black disk radius Rtot = tj mot /27r (dashed red curve) as functions of the center-of-mass energy W . 

The reason for this result can be seen in Fig. [S] where we plot the edge integrand normalized to Rtot- The resulting 
distributions have an approximately fixed shapes and areas as functions of FF, and migrate slowly to larger values 
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of b with increasing energy. We emphasize that these properties are determined up to 7 TeV by our ht to data on 
CTtot and (Tinei or CTeias- The fits to (Ttot extend to W ~ 70 TeV. The results shown in Fig.|S]at higher energies give 
our predictions based on the present eikonal model; the results are consistent with those of Block et al. [l[ which are 
independent of a detailed eikonal description of the scattering. 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

b, fm 

FIG. 9: Plots of the normalized edge integrands b{cosxR — /Rtot for, left to right, IT = 30 GeV (red), 500 GeV 

(green), 5000 GeV (brown), 5xl0'‘ GeV (blue), 10® GeV (purple), and 10® GeV (black). 

As seen in Fig. El tedge and i?tot cross in magnitude for FF ~ 2 TeV, with i?tot larger and increasing at higher 
energies. The eikonal amplitude for dtot is also beginning to saturate at 1 at small b in this region. The crossover 
point therefore gives a reasonable estimate of the energy at which the scattering amplitude begins to show aspects 
of asymptotic behavior, with the edge region becoming less important than the central region. We note that the 
crossover region is where the gluon contributions to the cross section become dominant as is evident in Fig. El 


D. Origin of the edge 

The two key results discussed above are: (1), the logarithmic growth with energy of the effective radius of the 
strongly absorptive black disk region in the scattering amplitude and the resulting In^ s growth of ctot; and, (2), the 
existence of an edge region with constant width tedge- 

There are several possible explanations for the origin and constancy of the edge in the pp and pp scattering 
amplitudes. A classic explanation would attribute the edge to pionic fluctuations around the proton or antiproton, 
with the pionic fluctuations then interacting in the collision. This would give an edge region on the scale of Ri 1.4 

fm, but does explain how the overall radial scale of the scattering amplitude would increase, or how strongly the 
fluctuations would couple to the the expanded proton in this picture. 

The most likely combined explanation for both effects in our view is connected to the rapid increase in the strength 
of the gluon helds in the interacting hadrons with increasing energy, and the resulting saturation of the exponentially 
bounded gluon confinement volumes of the two particles. This leads to a growing likelihood of interaction in a pp or 
pp collision and gluon dominance of the scattering, with the radius of the r^ion of strong absorption growing as Ins. 
Specific models based on parton collisions @ and Reggeon field theory [1,0] have this character. 

This picture should hold in any hadronic scattering at energies sufficiently high that the quark effects prominent 
at low energies can be neglected. Since the gluon confinement volume is presumably fixed, all hadronic cross sections 
should then have a universal constant x In^ s behavior at high energies, with a common value of the multiplicative 
constant. This behavior is consistent with the observed behavior of all cross sections which are known at high energies 

13 . 

The scale of the gluon conhnement volume is set in our model by 1/pgg. The increasing strength of the gluon fields 
is represented in QCD-based minijet-type models by the growth of the gluon distribution functions fg with decreasing 
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Bjorken x values or increasing energy of the gluon-gluon collision. This leads to stronger gluon-gluon scattering in 
the collision, typically increasing as a power of s, corresponding growth of xi, and gluon dominance of high-energy 
scattering. 

Our eikonal model and that of Fagundes et al. have this general character. The power-law growth of Ggg with 
s, combined with the exponential cutoff in the overlap function A{b,^gg) at large impact parameters, is sufficient to 
ensure an asymptotic In s growth o f Rtr,t a nd In^ s growth of the cross sections at very high energies independent of 
any appeal to the Froissart bound [1^, on the cross sections. (3^ The dominance of the gluons at high energies 
is evident in Fig.jSJ In the black disk limit, Ueiasi crinei —>■ O'tot/S, the leading In^ s terms in the cross sections cancel 
in Eq. (fT71) . The asymptotic parametrizations in Eq. (1141) and the equivalents for Ueias and (Ti„ei combined with the 
logarithmic growth of i?tot then indicate that tedge should be constant, or nearly so, at high energies as is observed. 

The width of the edge is also related to 1/fJ.gg in our model: the factor rj in the edge integrand 77(1 — 77 ), Eq. (flTl) 
increases from 0 to 1 over an interval in b proportional to l/fJ-gg centered around i?tot, while the factor (1 — 77 ) falls 
from 1 to 0 over a similar interval. This results in the edge integrands shown in Fig. [5] with the width of the peaks, 
hence tedge, proportional to \/^gg. 

A different, but potentially related, mechanism was proposal by Rosner who described the edge in terms of 
the breaking of flux strings connecting quarks or 33 gluon conhgurations in the two hadrons. His estimate of the edge 
width, based on the energy needed to break such QCD strings in other processes, is of the right general size. We 
attribute the scattering mainly to interactions of the gluons or gluon helds in the overlapping nucleons in the collision. 
Some components of those fields may be mixed between the nucleons in the scattering, leading to the stretching of 
flux tubes between the nucleons as they separate and extending the effective range of the interaction through string 
breaking as Rosner proposed. The growth in radius of the main interaction region, however, arises from the increasing 
saturation of the gluon confinement volume as sketched above. 

Since gluon dominance should appear at high energies in all hadron-hadron scattering, we expect all hadron-hadron 
total cross sections to approach a universal In^ s growth with a common coefficient at very high energies as noted 
above. All hadron scattering amplitudes should also have an edge region with an approximately constant width 
proportional to Xj[igg\ the leading terms in the cross section difference in Eq. (1211) and its analog for other hadrons 
cancel in the difference given their universal behavior at high energies, and the subleading terms are presumably 
logarithmic in s as is the factor in the denominator. 

The asymptotic mass scale and the corresponding behavior of the eikonal function at large impact parameters are 
clearly of considerable interest. The model we have used assumes that x( 6 , s) can be written as a sum of terms in 
which the energy dependence factors out of the overlap functions A(b, A) where the latter, exploiting ideas originally 
formulated by Wu and Yang [ 2 ^, are given as convolutions of density distributions similar to those associated with 
the proton charge and magnetic moment form factors. The resulting overlap functions are small at large impact 
parameters, so the integrand for the scattering amplitude /(s,t) is proportional to <Jgg[s) x fj,gg{^ggb)^K^^figgb) at at 
large b for energies where gluon scattering is dominant. 

This is inconsistent on the surface with the result expected from the dispersion relation in t for /(s, t), schematically 



a{s,t') 
t' -t ' 


The partial wave amplitude for angular momentum j is just 


fjis) = 



( 22 ) 


(23) 


Using the standard approximation Qj (z) m Kq ^\/2j^(z — 1) j , very good for j large and (z— 1) small, and introducing 
the impact parameter b = j jp, we get the impact parameter representation 


f] (s) —^ ^ ^ (bVt''^ . (24) 

The functions x^K^[x) and Kq{x) behave quite differently for x large, and Eq. (l24l) involves an integral over t' while 
the gg eikonal function involves only the fixed scale pgg. There is consequently some uncertainty as to how well the 
asymptotic behavior of the scattering amplitude is described at large b in the present model. Seen a different way, 
the weight function a(s, t') for the gg term is proportional to 5{t' — Pgg)/{t' — t)^, the result obtained from a product 
of dipole form factors consistent with the proton electric form factor. 

This remains an interesting problem which deserves further study. 
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E. The edge and diffraction dissociation 

An interesting connection between the edge and diffraction dissociation follows from an old analysis of the latter 
by Pumplin Q, who used an argument based on unitarity and the properties of scattering eigenstates to show that 
the 6-dependent partial cross sections for the dissociation of an incident particle on a nucleus were bounded above 
by (l/2)crtot(6) — creias(6) where the cross sections refer to particle-nucleus scattering. In the present case of strongly 
absorptive pp or pp scattering, this argument leads to an upper bound on the single-particle dissociation cross section 
for either incident particle, 


^SD ^ (^tot 2fJelas)/2 (2^) 

nOC 

= 2tt dbbr]{l — rj), (26) 

Jo 

or a total dissociation cross section ddiss < crtot ~ 2creias when both possibilities are included. The partial cross 
sections — the sum of squares of the dissociation amplitudes — are similarly bounded by the integrand in Eq. (1261) , 
crsoib) < 77(1 - 77 ). 

The expression in Eq. (E51) is just the edge integral, so the Pumplin bound relates aso to the area associated with 
the rim of width ledge in the pp or pp scattering amplitude. Since ledge is essentially constant as seen in Fig. |S] and the 
edge integrand is centered on i?tot = \/ o'tot/27r, the bound on <730 grows proportionally to ^/d^ and will increase 
logarithmically at large s. 


o’SD < constant X In s, (27) 

s—^oo 

a result which follows from the established In^ s growth of CTtot for s —> 00 . 

Because Eq. (EH) only gives an upper bound on the growth of asD with energy, and there is no comparable energy- 
dependent lower bound, it is not clear that the real dissociation cross section will actually grow at this rate; however, 
that seems likely given the growth of other hadronic cross sections. The possibility of this behavior was noted in 
[ 2 ^ for a specific model of diffraction dissociation, but some earlier models predicted a decreasing dissociation cross 
section. 

Given the observed constancy of the edge width and the bound in Eq. (1251) . we find that the ratio of the dissociation 
cross section to the total cross section must decrease at least logarithmically at high energies, o’sD/o’tot < constant x 
1 /.Rtot oc 1 / In s —>■ 0 for s —>■ 00 . 

As emphasized in [^ . many studies of particular mechanisms for diffractive dissociation neglect aborptive effects. 
These are clearly crucial in Eq. (1261) : 77 is very small at high energies out to impact parameters 6 near Any 

reasonable model of diffractive dissociation must take this into account. 

We have calculated the bound on the single particle dissociation cross section using the expression in Eq. (l26l) and 
the eikonal factor 77 ( 6 , s) found in our ht to the combined pp and pp data. We show the result and the CDF [s^ 
measurements of diffractive dissociation in this process in Fig. 1101 

As expected, the measured cross section for diffractive dissociation lies considerably below the calculated inclusive 
value: the upper limit on asD is only reached under special conditions Furthermore, most experiments concentrate 
on the differential cross section dJasn/dtdM^ in order to test specific models. The kinematic regions in which this 
can be measured and distinguished experimentally from other inelastic processes, are quite limited. It is typical to 
require, for example, a very large ratio of the final center-of-mass momentum p' of the surviving particle to its initial 
momentum p, p'/p > 0.85 in the CDF experiments [s^, and further conditions on the detectability and mass of 
the dissociated system. There appear, in fact, to be no universally accepted experimental criteria for extracting this 
cross section, with the results typically depending on how the distributions in the momentum transfer t and Mj- 
are modeled. The CDF results integrated over the allowed regions give cross sections which do increase with energy 
roughly as predicted by the bound as shown in Fig.fTOl 

Finally, to get an idea of the expected dependence of the scattering on the momentum transfer, we have calculated 
the analog of the elastic scattering amplitude /(s, t) , Eq. ([T|) , for the exclusive process p-\-p ^ p + X using the bound 
77(1 — 77 ) on the 6 -dependent amplitude noted above, 

pOO 

/sd(s,(?)=/ dbbr]{l-■r])Jo{qb). (28) 

Jo 

Here = 2pp'{l — cos 9) = —t — — 777^)/2 + ■ ■ ■, W ^ where 9 is angle through which the the surviving 

particle is scattered and Mx is the mass of the system X [^. We do not specify the dependence of the dissociation 
process on Mx, but think of /sd as giving the characteristic 6 dependence of the mass-dependent amplitudes averaged 
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FIG. 10: The curve gives the upper bound on the inclusive cross section for single-particle diffractive dissociation p+p p + X 
calculated using ctsd < (ftot — 2geias)/2 and our eikonal fit to the pp and pp cross section data from 5 GeV to 70 TeV. The 
data shown for a so are from GDF 


over masses. The slope Bsd of the corresponding cross section at 0 = 0 is defined for purely absorptive scattering 
{XR = 0) as 

nOO pCO 

Bsd = dln\fsD\‘^/dq‘^ = (1/2) dbb^r]{l - ?]) / / dbb7]{l-r]). (29) 

^0 ' do 

The results we obtain for |/sd|^ are shown at several energies are in Fig. ITT] 

We find that the cross section corresponding to the smooth peripheral edge distribution in Eq. (1251) develops diffrac¬ 
tion zeros at a given energy W at smaller values of than the elastic cross section and that the forward slope 
parameter is larger, Bsd > ^eias- For example, at W = 1000 GeV, Bsd = 23.9 GeV“^ while Belas = 16.0 GeV“^. At 
100 GeV, Bsd = 18.0 GeV“^ while Beias = 13.3 GeV“^. We emphasize that these results assume that the specific 
mechanism in question for p+p ^ p + X is represented in impact parameter space by a distribution oc rj {1 — rj) which 
covers most of the region allowed by Eq. (1261) at the given energy. 

In the opposite extreme in which the mechanism for the dissociative production of a particular system X is repre¬ 
sented by a narrow distribution around an impact parameter bo in the allowed region, /sd cx; Jo{qbo), and Bsd = bo/2. 
Eor bo near the peak of the edge distribution, bo ~ i?tot = \/ o’tot/27r, the slope parameter is somewhat smaller than 
Belas, 13.9 GeV“^ instead of 16 GeV“^ atW = 1000 GeV, and 9.3 GeV“^ versus 13.3 GeV“^ at IF = 100 GeV. This 
picture is general because of the compact nature of the dissociation distribution in Eq. (1261) . Any model for, or mea¬ 
surement of, diffraction dissociation must give a slope parameter in the range spanned by these limiting cases, most 
likely slightly larger than B//^^/2 because the extra factor of b^ in the numerator of Eq. (I29|) weights that distribution 
toward larger impact parameters than that in Eq. (1261) which is centered at ~ Rtot ■ 

Finally, we emphasize that it would be of considerable interest to measure the total single-dissociation cross section 
CTSD, for example for p + p ^ p + X, without strong restrictions on p' and Mx and only the requirement that the final 
state contain an isolated p near the forward direction at a fairly low momentum transfer q, opposite a multi particle 
system. This would determine how close the bound on asD is to saturation. 


F. Survival of rapidity gaps 

The search for new physics in pp or pp collisions can be simplified when the new process occurs in a large rapidity 
gap so is not accompanied by unrelated secondary particles in that region. An example discussed by Block and Halzen 
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FIG. 11: The model differential dissociation cross sections calculated as the squares of the amplitudes /sd in Eq. (1281) for 
W — 546 GeV (black dashed curve), 1800 GeV (solid blue curve), and 7000 GeV (dot-dashed red curve). The results illustrate 
the and W dependence expected for dissociation ampitudes which saturate the edge distribution in impact parameter space, 
but are not predictions for the actual Mx-dependent cross sections. 


would be Higgs boson production through W boson fusion, WW —1 H, where the W s are emitted by quarks in 
the colliding hadrons. 

We take the inclusive differential cross section for this process in impact parameter space as 

= crww^nMb, ^J-qq) (30) 

where A{b, ^qq) describes the spatial overlap of the quarks distributions as defined in Eqs. (IA3I) and (IA4|) . 

A{b, fiqq) is normalized so that integration over (Ph with no further input would just give aww^H as calculated in 
the parton model. However, further inelastic processes can occur in the hadronic collision giving secondary particles 
other than those associated with the remnants of the incident particles, and eliminating the rapidity gap. The 
probability that no such inelastic process occurs is , and the cross section including this survival probability is 
therefore 


= o-ww^HA{b, ^iqq)e (- 32 ) 

The factor A(6,is just the differential survival partiality for the gap. The construction generalizes to 
other processes. 

Defining this following [3l| as d{\S\^)/(Ph, the total survival probability for the gap is 

(|5|")= J d^hAib,^lqq)e-^^^^^’^'>. (32) 

These survival probabilities were calculated in [3l| for the eikonal model discussed there. We have recalculated the 
survival probabilities using the eikonal model developed here. The results, given in Table HI are very similar. The 
same calculation is easily done for gluon-initiated processes. 


V. SUMMARY 


In this paper, we presented the results of a detailed analysis of the current data on pp and pp scattering in the 
eikonal formalism, parametrizing the eikonal function in a form suggested, but not restricted, by the structure found in 
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TABLE I: The gap survival probabilities (|5'|^) in percent for pp and pp collisions as functions of the center-of-mass energy W. 


W, GeV 

pp {%) 

PP (%) 

63 

38.7 ±0.6 

38.4 ± 0.6 

546 

28.6 ±0.5 

28.6 ± 0.5 

630 

27.8 ±0.5 

27.8 ± 0.5 

1,800 

22.2 ±0.5 

22.2 ±0.5 

14.000 

13.1 ±0.3 

13.1 ±0.3 

40,000 

9.8 ±0.2 

9.8 ±0.2 


minijet models for the scattering. The fit to the combined data is excellent. Our results for the total and elastic cross 
sections and p values agree very well with the earlier fits of Block and Halzen based on analytic amplitudes 

with a In^ s growth at high energies. 

We showed that our model, which includes a gluonic contribution to the eikonal function that grows as a power of 
s = leads naturally to cross sections, p values, and slope parameters B which can be described very accurately 
at high energies by quadratic expressions in Ins. We can therefore extend the earlier Block-Halzen analysis of the 
energy dependence of the total cross sections and p values to include the elastic and inelastic scattering cross sections 
and the logarithmic slope B of the forward elastic scattering cross section, all measured quantities. 

Our detailed model allowed us to analyze the impact-parameter structure of the various scattering amplitudes in 
detail, including the relative importance of various contributions to the amplitudes, and the approach to asymptotic 
behavior at high energies where gluonic processes become dominant. We commented on uncertainties in the asymptotic 
behavior of the eikonal function which affect the asymptotic behavior of the cross sections and merit further study. 

We used our model to examine the structure of the edge of the pp and pp scattering amplitudes recently identified 
by Block et al. [l[ in some detail. The width of this edge region is nearly energy independent at ~ 1 fm, a property 
clearly evident in our results. We commented on some possible ^namical origins for the edge. 

We also used the model to investigate the Pumplin bound Q on the cross section for single particle diffractive 
dissociation which is given directly in terms of the edge cross section. The constancy of fedge and the In^ s growth of 
CTtot at high energies show that the bound — and possibly asD — increase only as Ins, while the ratio cfsd/ crtot must 
decrease at least as 1/lns. 

Finally, we used the model to update earlier results on the survival probability of large rapidity gaps in pp and 
pp scattering, a matter of interest in the search for rare processes in the scattering. 
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Appendix A: The modified Aspen model 

The model used here in fitting the pp and pp cross sections is a modification of the “Aspen model” of Block et al. 
[i,®- That model was based the structure of the eikonal function found in QCD minijet models for the scattering 
in which the interactions between hadrons are described in terms of the interactions of their constituent quarks and 
glu ons with allowance for “soft” interactions at low momentum transfers. While we will follow the notation used in 
[lOl| . the identification of the terms made there as describing quark-quark {qq), quark-gluon (qg), or gluon-gluon {gg) 
interactions becomes blurred in the general setting, especially for the qq and qg terms. 

We will write the eikonal functions in terms of crossing-even and crossing-odd components as 

XEib,W) = i[(jgg{w)A{b,pqq) + aqg{w)A{b,fj.qg)+agg{w)A{b,fj.gg)], 

Xo{b,W) = -J:ggC5['^y~^"^ A{b,Podd), 


(Al) 

(A2) 
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where one needs to make the replacement w ^ We in the final results to obtain the correct asymptotic phase 
required by analyticity and crossing symmetry i[ni- Here we will simply write the functions on the right-hand 
sides of Eq. (lAlIl and Eq. (IA2I) as functions of rc, with the replacement to be made in the final results. The constant 
Egg = 97 ra^/mg sets the scale in Eq. (IA2I) and later equations. 

The overlap factors A{b, A) in these expressions are defined in terms of the relevant distributions in the proton by 

Aij = j d?hpi(h')pj{h-h'). (A3) 

Assuming that the distributions pi have approximately the same form as that determined from the proton electric 
form factor, the overlap functions become 

^2 rco 

A{b,X) = —{XbfK3{\b), d^hA{h,X) = l (A4) 

for appropriate choices of the A parameters. 

The gluon-gluon term in Eq. m, dominant at very high energies, was parametrized in 0 using a very simplified 
description of gg scattering in low-order QCD. The result was an expression which involved a leading power of s/itiq, 
logarithms of that quantity, and a constant term, plus terms involving inverse powers of sjrn^. Given the uncertainties 
in the model, including a rather arbitrary choice of the leading power, we will simply parametrize Cgg directly in terms 
of a power and leading logarithm in s/mg with an additive constant chosen so that the gg term gives a negligible 
contribution to the eikonal function at low energies as in low-order QCD. 

The gg and qq terms have a less singular structure in QCD, and model results derived using scaling parametrization 
of the quark structure function /g do not separate cleanly from the expected contributions from soft processes or 
the gg terms. We simply follow the parametrizations used in 0 allowing, however, the powers in the Regge-like 
low-energy terms to vary from the ^/^/s = XjW behavior assumed there. 

The “cross sections” Uij in Eq. (lAll) are then 

Oqqiw) = Egg [Co-b , (A5) 

aqg(w) = EggC 2 ln(wVmg), (A 6 ) 

agg(tc) = Egg {0.0713+ C3ln(W/W'o) + C4[(tE/mo)'^-(Wo/™o)^]} (A7) 

The expression in Eq. (IA7I) gives an excellent fit to the more complicated and restricted form for Ugg derived in 0 
at higher energies where this term is important, and contributes less than 1 % of the total eikonal function at b = 0 
at the low-energy matching point W = Wq = 4 GeV in agreement with the results there. 

Our fit to the pp and pp scattering cross sections, the ratios p of the real to the imaginary part of the forward 
elastic scattering amplitudes, and the logarithmic derivatives B of the forward differential cross sections dcr/dt used 
the six coefficients C^Ci, C 2 , C 3 , C 4 , C 5 and the parameters ai, a 2 , and [3. The remaining parameters /Xgg, pgg and 
mo were fixed as in [lOj with the energy scale mo = 0.6 GeV, and the ps chosen by hand in the range determined by 
the proton charge form factor, pgg = 0.705 GeV and pgg = 0.89 GeV. We did not vary these parameters in making 
the fit. We note also that the overall factor Egg = 97 rag/mg which appears in the cross sections in Eq. (IA2I) and Eqs. 
(IA5IIA5I) can be absorbed into the coefficients Cp it was separated out in ii to provide a connection with minijet 
models for the eikonal function y where such factors appear naturally. 

A summary of the parameters with the results of the fit is given in Table |TT1 
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